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flavor-dependent angular distributions 
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The usual description of self-induced flavor conversions for neutrinos (v's) in supernovae is based 
on the simplified assumption that all the p's of the different species are emitted "half-isotropically" 
by a common neutrinosphere, in analogy to a blackbody emission. However, realistic supernova 
simulations show that f angular distributions at decoupling are far from being half-isotropic and, 
above all, are flavor-dependent. We show that flavor-dependent angular distributions may lead to 
crossing points in the angular spectra of different v species (where _F^^ = F^^ and Fj^^ = F^^ ) 
around which a new multi-angle instability can develop. To characterize this effect, we carry out 
a linearized flavor stability analysis for different SN neutrino angular distributions. We find that 
this instability can shift the onset of the flavor conversions toward low-radii and produce a smearing 
of the splitting features found with trivial i/ emission models. As a result the spectral differences 
among i/'s of different flavors could be strongly reduced. 

PACS numbers: 14.60. Pq, 97.60.Bw 



Introduction. — The neutrino flux emitted from a core- 
collapse supernova (SN) represents a powerful tool to get 
valuable information about the mixing parameters and 
the dynamics of the exploding stellar core. In this con- 
text, renewed attention is beiiig paid to collective fea- 
tures of flavor transformations [l|, |2| induced by ly-i/ self- 
interactions d, i] in the deepest SN regions, near the 
neutrino-sphere (see for a recent review). In partic- 
ular, it has been claimed that the self-induced "spectral 
splits" would be observable in the v burst from the next 
galactic SN, allowing to get crucial information about the 
unknown v mass ordering (see, e.g., @). 

The development of these self-induced effects crucially 
depends on the inner boundary conditions fixed for the 
flavor evolution, e.g. on the v initial energy and angu- 
lar distributions. In this context, all the recent numer- 
ical simulations have been based on the so-called "bulb 
model" (see, e.g, [H, Q)- Assuming that the v trapping 
region at high matter density can be decoupled from the 
region at lower density where flavor conversions would 
start, v^s of different species are considered as emitted 
"half-isotropically" (i.e. with all outward-moving angu- 
lar modes equally occupied and all the backward-moving 
modes empty) by a common spherical "neutrinosphere," 
in analogy with a blackbody emission. However, this sim- 
plified toy-model may not capture important features of 
the SN v emission. In particular, the transition between 
the isotropic v emission in trapping regime (at a higher 
matter density) to the forward-peaked free streaming (at 
a lower matter density) docs not necessarily imply a ra- 
dius where the v angular distributions are half-isotropic 
(see, e.g., 0, 111)- Moreover, the different v species de- 
couple at different radii. Therefore, the physical last- 
scattering neutrinospheres for the different species would 
not coincide. Then, fixing a common neutrinosphere as 
inner boundary for the flavor evolution, one would real- 
istically find there different angular distributions for the 



different v species. In particular, since the Fe's and i^x's 
decouple at smaller radii with respect to i^e 's their distri- 
butions are more forward-peaked. 

The presence of non-trivial angular distributions was 
claimed in to produce a multi-angle instability in 
the self- induced flavor evolution of a toy model of v gas. 
However, in that seminal paper the author warned read- 
ers not to draw firm conclusions for the realistic SN 
case from his analysis performed with a small number 
of angular modes, before having explored this case with 
large-scale numerical simulations. Stimulated also by his 
intriguing hint, we decided to begin a systematic study 
of this issue. We find that as long as flavor universality 
is not broken in the v angular distributions, only minor 
quantitative effects appear in the evolution with respect 
to the half-isotropic case. Conversely, flavor-dependent 
angular distributions can dramatically affect the conver- 
sions of the dense SN i' gas. We relate this behavior 
to the presence of crossing points in the angular spectra 
of different v species (i.e., at emission angles for which 
F,y^ = F,y^ and = ^p^)- It is known that ly energy 
spectra presenting crossing points in the energy variable 
can develop instabilities around them This effect 

has been analytically understood performing a linearized 
stability analysis of the v equations of motion [l2|. Here, 
applying the same analysis, we find that an analogous ef- 
fect occurs when the angular spectra of different flavors 
present crossing points. As a consequence, this multi- 
angle instability can lead to a wash-out of the splitting 
features found assuming a half-isotropic emission. 

Energy and angular distributions. — In order to per- 
form our study of the multi-angle instability we have 
to fix our ly emission models. Our main goal would 
be to remove the bulb-model approximation of a half- 
isotropic neutrino angular distribution at the neutri- 
nosphere, equal for all the species. In order to introduce 

angular distributions, it is convenient to parameterize 
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FIG. 1: Upper panels: Multi-angle flavor evolution for //'s 
(left panel) and F's (right panel) . Energy spectra initially for 
Ve (black continuous curves) and Vx (light continuous curves) 
and after collective oscillations for i/e with half-isotropic an- 
gular distribution (dashed curves) and for an angular distri- 
bution with /3e = 1-0, Px = 1-5 (thick dotted curves). Lower 
panels: Difference of the energy-integrated angular spectra 
Fv^{u) — Fv^{u) for i/'s (left panel) and F's (right panel) at 
the neutrinosphere for the half-isotropic case /3e = = 0.0 
(dashed curves) and for the cases with fie — 1.0, Px = 1.5 
(dotted curves), j3e = 1.0, fix = 3.0 (continuous curves) and 
Pe = Px ~ 3.0 (dash-dotted curves). 



every angular mode in terms of its emission angle iJjj rel- 
ative to the radial direction of the neutrinosphere, that 
here we sehematically fix at i? = 10 km. A further sim- 
plification is obtained if one labels the different ang ular 
modes in terms of the variable u = sin^ "dn, as in jl2l. [l3|. 
Note that for an half-isotropic emission at the neutri- 
nosphere the ly angular distribution of the radial fluxes is 
a box spectrum in < m < 1 fl^ . 

We remark that energy and angular distributions of 
SN v's entering the ly number fluxes F^^{E,u), inte- 
grated over a sphere of radius r, are not independent 
of each other. However, here we schematically assume 
that the angular distributions are energy independent. 
Then, we can factorize the v flux of each flavor as 
F^^ {E, u) = N^^ X (f^^ {E) X C/^„ {u). The function (p^^ {E) 
is the normalized v energy spectrum (J dE(p^^{E) = 1), 
which we will parametrize for the different flavors with 
deformed Maxwell-Boltzmann distributions, as suggested 
by SN simulations [14|. The function U^^{u) is the nor- 
malized angular distribution {J duU^^{u)=l). We fix 
the neutrino average energies at {{E^^), (Ep^), {Ey^)) = 
(12, 15, 18) MeV. Concerning the possible v flux or- 
dering we consider a case representative of the cooling 
phase [3, i.e. N^^ : N^^ : N^^ = 1.13 : 1.00 : 1.33. 
According to this choice, in the upper panels of Figure 1 



are represented the initial fluxes for v^. (continuous black 
curves) and (continuous light curves) for v's (left pan- 
els) and D's (right panels) respectively. 

A systematic characterization of the v angular distri- 
butions at different radii in SNe is lacking in literature. 
Therefore, inspired by an inspection of some real an- 
gular distributions |15|, we use a simple toy model to 
capture the main deviations with respect to the half- 
isotropic bulb model, where U^^ = 1 for all the v species. 
In particular, we choose forward-peaked distributions 
Ui,^{u) cx (1 — u)'^°/^. For simplicity in the follow- 
ing we assume U,^^ = fJp^. In the lower panels of Fig- 
ure 1 we plot the energy-integrated spectral difference 
Fi,^ {u) — F,y^ (li) for r/'s (left panels) and the analogous one 
for 77's (right panels) for four representative (/3e, f3x) cases. 
Namely we compare the half-isotropic case (/3e = fix = 0, 
dashed curves) with three non-trivial angular distribu- 
tions: /3e = 1.0, Px = 1-5 (dotted curves), /3e = 1-0, 
/3x = 3.0 (continuous curves), and P,, = j3x = 3.0 (dash- 
dotted curves). One realizes that in the half-isotropic 
case /3e = 13 x = 0.0 and in the flavor universal case 
Pe = Px = 3.0, the differences of angular spectra do not 
present any crossing point in the angular variable, while 
in the other two cases, the spectra present a crossing 
point at finite < w < 1, where Fy^{u) = Fy^{u). We 
will show that crossing points in the energy-integrated 
angular spectra will produce a significant speed-up of the 
multi-angle instability in the v flavor conversions. 

Setup of the flavor evolution. — Our description of 
the non-linear v flavor conversions is based on a two- 
flavor oscillation scenario. For simplicity, we refer to the 
cases studied in [l^ in which three-flavor effects associ- 
ated with the solar sector are sub-leading. In this sit- 
uation, self-induced oscillations are driven by the atmo- 
spheric mass-square difference /S.rn^^^ = mg — m\ 2 — 
2.0 X 10^'^ eV^ and by a small (matter suppressed) in- 
medium mixing Om = 10""^ We refer to late-time 
cooling phase (t > 1 s) where the electron density rif. is 
smaller than the v one rii, and then has a subleading 
role on the development of the collective oscillations [l^ ■ 
Apart from the mixing suppression, we neglect the sub- 
leading matter effects in the flavor evolution. We will 
refer to the inverted mass hierarchy (Am^^j^^ < 0). 

Following [1^ . we write the equations of motion for 
the flux matrices ^e,u as function of the radial coordi- 
nate. The diagonal elements are the ordinary num- 
ber fluxes Fi,^{E,u). We normalize the flux matrices 
to the total Vg number flux Np^ at the neutrinosphere. 
Conventionally, we use negative E and negative num- 
ber fluxes for D^s. The off-diagonal elements, which are 
initially zero, carry a phase information due to flavor 
mixing. Then, the equations of motion read idr^E,u = 
[He,u,^e,u] with the i^-iy Hamiltonian [sl. R [T2j 
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The factor proportional to the neutrino velocity Vu,r 
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(1 — ui?^/r^)^''^ [T^ implies "multi-angle" effects for neu- 
trinos moving on different trajectories [H, d, |3l- In order 
to properly simulate numerically this effect one needs to 
follow a large number [0(10'^)] of interacting v modes. 

Stability condition. — In order to achieve a deeper un- 
derstanding of the multi-angle instability, triggered by 
non-trivial angular distributions, we find particularly 
useful to apply to our problem the linearized stability 
analysis recently worked out in [l^l- This method al- 
lows us to determine the onset of the flavor conversions, 
seeking an exponentially growing solution in the eigen- 
value equations, associated with the linearized equations 
of motion for the v ensemble. 

In order to perform the stability analysis we closely 
follow the prescriptions presented in |12| and summa- 
rized in the following. At first we switch to the fre- 
quency variable a; = Am^^jjj/2£' and we introduce the 
neutrino flux difference distributions g^^^u = 5(w,u) = 
\dui / dE\[Fy^ {E{uj), u) — F^^ {E{u:),u)\, normalized to the 
total flux at the neutrinosphere. Note that g^j^u is 
defined also for negative oj, where it represents the dif- 
ference of fluxes in the antineutrino sector in th e op posite 
ordering. Then, we write the flux matrices as jl2| 



Tr$„ 



9lo,u 



5, 



UJ,U 



(2) 



where Tr is conserved and then irrelevant for the fla- 
vor conversions, and the initial conditions for the "swap- 
ping matrix" in the second term on the right-hand side 
are Su:,u = 1 and S^^u = 0. Self-induced flavor transi- 
tions start when the off-diagonal term Suj,u exponentially 
grows. In the small-amplitude limit l^t^,?!! ^ 1, and at 
far distances from the neutrinosphere r ^ R, the lin- 
earized evolution equations for Su,.u lead to an eigenvalue 
equation [l^l 



{oj+ue^~VL)Qt^^u ^ 1^ I du'duj' {u+u')g, 



(3) 



for the eigenvector Quj,m obtained writing the solution 
of the hnearized equation S^^^u = Qu,ue~^^^^. Here we 
introduced the complex frequency — j + in and the 
parameter e = J du du gi^^u quantifying the "asymme- 
try" of the neutrino spectrum, normalized to the total 
Ve number flux. For our specific choice of v spectra, it 
results e = 0.13. The I'-i' interaction strength is given by 



2r2 



(4) 



For our SN model, it results ^{R) = 4.7 x 10"'' km^^ A 
solution of Eq. ([3]) with k > would indicate an expo- 
nential increasing Quj,u, i-C. an instability. When an in- 
stability occurs, for a given angular mode uq the function 



is a Lorentzian [111, centered around a resonance 



frequency uj ~ j ~ uefj,, and with a width k. 

Flavor conversions — In order to illustrate the effect 
of the angular distributions, in Fig. 1 (upper panels) we 
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FIG. 2; Upper panel: Radial evolution of the integrated z- 
component of the polarization vector Pz for i^'s for the half- 
isotropic case (dashed curve) and for angular distributions 
with /3e = 1.0, 13^ = 1.5 (dotted curve), pe = 1.0, 13^ = 3.0 
(continuous curve), and — I3x ~ 3.0 (dot-dashed curve). 
Lower-panel: Radial evolution of the eigenvalue ft for the same 
cases shown in the upper panel. 



compare the oscillated electron (anti)neutrino fluxes (at 
r = 300 km) in the half- isotropic case (/3e = /3x = 0', 
dashed curve) with a non-trivial angular distribution 
with /3e = 1.0, l3oo = 1.5 (thick dotted curve). The ef- 
fect on the flnal spectra is dramatic: The spectral splits 
observable in the flnal i^e spectra in the half-isotropic case 
are smeared-out and the spectral swaps are not complete. 

In the upper panel of Fig. [2] we show the radial evolu- 
tion of the z-component of the integrated value of the 
polarization vector P^, that is related to the flavor con- 
tent of the ensemble for four different angular distribu- 
tions of Fig. [T] Namely, we compare the half-isotropic 
case (/3e = 0, fix = 0, dashed curve) with three non- 
trivial angular distributions: /3e = 1-0, f3x = 1-5 (dot- 
ted curve), /3e = 1-0, I3x ~ 3.0 (continuous curve), and 
Pe = Px ~ 3.0 (dot-dashed curve). First of all, note 
that with the /3e — fix = 3.0 flavor-blind distributions 
the behavior does not present major changes with re- 
spect to the naive case, except from a delay in the start- 
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ing of flavor conversions at r ~ 90 km. In the case of 
Pf. = 1.0, = 1-5, the onset of the flavor conversions 
(at r ~ 70 km) is very close to what expected in the 
half-isotropic case (r ~ 75 km). However, the further 
flavor evolution shows dramatic differences. Indeed, with 
the non-trivial angular distribution the final value of the 
integrated v polarization vector would be ~ indi- 
cating that the spectral differences between the final Vf, 
and Vx spectra are strongly reduced with respect to the 
half-isotropic case. Choosing larger differences in the an- 
gular distributions of different fiavors, as in the case with 
Pe = 1.0, Px = 3.0, not only the final value of the Pz will 
change, but also flavor conversions would start much ear- 
lier (r ~ 55 km) than in the half-isotropic case. 

In the bottom panel of Fig. [2] we plot the radial evo- 
lution of the eigenvalue k, coming from the solution of 
Eq. Q for the same cases shown in the upper panel. In 
the half-isotropic case Pe = Px = ^ (dashed curve) the 
K function presents a hump peaked around r ~ 100 km 
and connected with a long tail around r ~ 75 km. This 
tail indicates that the system is in principle always unsta- 
ble. However, one can verify that in the tail the unstable 
v frequency modes are in the infrared region, where the 
spectrum is strongly suppressed. Therefore, they have no 
impact for the flavor conversions. The onset of the fla- 
vor conversions is then given by the connection between 
the hump and the tail at r = 75 km. This result is in 
agreement with what shown in the numerical calculation 
of Pz, shown in the upper panel (dashed curve). 

When considering non-trivial angular distributions, 
the long tail in the k function in the half-isotropic case 
now disappears. However, the presence of a non-isotropic 
angular distribution is not enough to produce an en- 
hancement in the value of k. Conversely, in the flavor 
blind case P^. = Px = 3.0 (dot-dashed curve) the insta- 
bility is suppressed with respect to what has been seen 
in the half-isotropic case. In this case, flavor conver- 
sions start at r = 80 km. A signiflcant enhancement of 
the multi-angle instability occurs when v angular spec- 
tra exhibit crossing points in u, as can be seen in the 
case with /3e = 1.0, Px = 1-5 (dotted curve) and even 
more in the one with P^ = 1.0, Px = 3.0 (continuous 
curve). In both cases the peak in k increases with re- 
spect to the half-isotropic case and also the hump broad- 



ens toward smaller r. Since k reaches a higher peak 
value, the width of the Lorentzian around an unstable 
frequency mode for a given angle u would be broad, im- 
plying a speed-up in the transitions. In particular, in 
the case of Pe = 1.0, /S^, = 1.5 flavor conversions start at 
r ~ 70 km, while in the case with /3e = \Xs,px = 3.0 
around r ~ 55 km, in agreement with the numerical re- 
sults shown in the upper panel of Fig. [2] 

Conclusions — We have described a new instability in 
self-induced flavor conversions for SN i^'s, associated with 
flavor-dependent angular distributions. These can lead 
to crossing points among the different spectra that can 
produce new flavor conversion effects, absent with an 
half-isotropic f emission. We checked that this effect 
would develop in both the mass hierarchies and would 
be particularly relevant during the SN cooling phase, 
where the differences among the different flavors are rel- 
atively small, and self-induced effect can develop with- 
out any matter hindrance. The effect on the oscillated 
V spectra can be dramatic. Namely, it would produce a 
smearing of the splitting features widely discussed in the 
half-isotropic case, and resulting v fluxes with less sig- 
nificant spectral differences. This tendency toward spec- 
tral equalization would challenge the detection of further 
oscillation signatures, like the ones associated with the 
Earth crossing of SN z^'s (see, e.g., [13 )• Moreover, we 
found that also onset of the flavor conversions can be 
signiflcantly pushed at low-radii, challenging the multi- 
angle suppression found in the half-isotropic case [20| . 
Possible flavor conversions at low-radii would have an 
interesting impact on the r-process nucleosynthesis in 
SNe [111 . Furthermore, if i/ oscillations develop too close 
to the neutrinosphere, they would invalidate the v trans- 
port paradigm in SNe that ignores v conversions. Given 
the importance of the effects discussed here, a future task 
would be to perform the flavor evolution and stability 
analysis taking into account angular distributions as re- 
alistic as possible, obtained directly from SN simulations. 
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